£ N 8 1 - 1 4 6 9 9 ^ 

... 

Grid Generation for Time Dependent Problems: Criteria and Methods 

Marsha Berger, William Gropp and Joseph Oliger 
Department of Computer Science 
Stanford University 

Abstract 

We consider the problem of generating local mesh refinements when 
solving time dependent partial differential equations. We first discuss the 
problem of creating an appropriate grid, given a mesh function h defined over 
the spatial domain. A data structure which permits efficient use of the 
resulting grid is described. Secondly, we show that a good choice for h is 
an estimate of the local truncation error, and we discuss several ways to 
estimate it. We conclude by comparing the efficiency and implementation 
problems of these error estimates. 


WHAT ADAPTIVE MESH GENERATION FOR TIME DEPENDENT PDE'S 

OBJECTIVES REDUCE # MESH PTS 
MINIMIZE OVERHEAD 

TRADEOFF: EXTRA PTS. VS. EXTRA LOGIC 

REQUIREMENTS ion awia sovj owasoaad 

• MARCHING ALGORITHMS WILL BE USED 

• COMPUTING TRANS IENT SOLN BY FINITE DIFF. 

• TIMESTEP SMALLER ON FINER GRIDS , MESH RATIO CONSTANT 

• GRIDS MUST CHANGE WITH TIME 

• COARSEST GRID DOES NOT CHANGE WITH TIME 
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DESCRIPTION OF GRIDS 


LOCALLY UNIFORM 


RECTANGLES OF ARBITRARY ORIENTATION. EXTENSIONS TO 
CURVILINEAR GRIDS FITS INTO SAME FRAMEWORK 

SUPPOSE BASE GRID G n = V G A . FORM HIERARCHY OF 

u j u . J 

NESTED GRIDS WHERE EACH REFINED GRID IS WHOLLY 
CONTAINED IN A SINGLE COARSER GRID 



• REFINED GRIDS CAN BE CONSTRUCTED AUTOMATICALLY 
AT t = 0 FROM INITIAL DATA. 


HOW GRIDS ARE FORMED 

GIVEN A "MESH FUNCTION" h(s.y) USED TO DETERMINE 
WHERE TO PLACE REFINED GRIDS. 

FLAG GRID PTS. WHERE h(x,y) >e. 

| — CLUSTER 
• ORIENTATION 
• GOOD FIT? 



CLUSTERING 



• NEAREST NEIGHBOR 

d(PT. , CLUSTER) < dmax => 

PT. e CLUSTER 

• SPANNING TREES 

CONNECT ALL PTS. ACCORDING TO SOME CRITERIA. 
BREAK LONGEST LINKS, 

• MINIMAL SPANNING TREES 

• MINIMUM DIAMETER TREES 


ORIENTATION 

• FIT ELLIPSE TO FLAGGED PTS. OF A CLUSTER USING 1ST AND 
2ND MOMENTS. 

• USE MAJOR AND MINOR AXES OF THE ELLIPSE TO GET RECTANGLE 
ORIENTATION. 

(REF: D. GENNERY, "OBJECT DETECTION AND MEASUREMENT USING 
STEREO VISION") PROCS. IJCAI. 1979, pp 320 - 327 

• FIT MIN. BOX TO INCLUDE FLAGGED PTS. + SMALL BUFFER ZONE 
FOR SAFETY. 



GOODNESS OF FIT 

• RATIO OF FLAGGED TO UNFLAGGED PTS. 

• IF TOO LOW, RECLUSTER AND REFIT. 




KEEPING TRACK OF GRIDS 

NESTING SUGGESTS USE OF TREE STRUCTURE 


(REF: KNVTH, "ART OF COMPUTER PROGRAMMING". VOL 1) 



1) GRID LOCATION 

2) SPATIAL AND TEMPORAL STEP SIZES 

3) SIZE OF GRID 

4) 3 TREE LINKS 

5) PTR. TO INTERSECTING GRIDS 

6) MAIN STORAGE AREA PTR. 


POINTS TO NOTE 


1) EASY TO HANDLE FAIRLY GENERAL REGIONS. 

ALL THE WORK IN SETTING UP THE PROBLEM 
IS IN SPECIFYING THE LOCATION OF THE COARSE 
GRID AND ITS CONSTITUENT RECTANGLES. THE 
REST IS AUTOMATIC. 




2) EASY TO USE DIFFERENT METHODS ON DIFFERENT 
GRIDS. 



WHAT IS h(x,y)? 


WOULD LIKE TO EQU I DISTRIBUTE THE GLOBAL ERROR. 

ID LINEAR THEORY SAYS IF CONTROL 

(1) INITIAL ERROR 

(2) BOUNDARY ERROR 

(3) LOCAL TRUNCATION ERROR 

AND METHOD IS STABLE FOR IBVP THEN THE METHOD CONVERGES. 

(1) AND (2) CONTROLLED BY STD. MEANS 
(3) CONTROLLED BY REFINING MESHES 

USE LOCAL TRUNCATION ERROR FOR h(x,y). 


REQUIREMENTS FOR LOCAL ERROR ESTIMATOR 

• ACCURATELY MIMIC ERROR BEHAVIOR 

• REASONABLY ACCURATE ESTIMATE - NOT NEC. A BOUND 

• CHEAP TO COMPUTE 

FLEXIBLE - EASY TO SWITCH INTEGRATORS 


• THE FEWER TIME LEVELS THE BETTER. 


POSSIBLE ESTIMATORS 


DIRECT ESTIMATION OF TRUNC. ERROR 
• FIND LEADING TERM 
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( e .g. - V t 
ESTIMATE BY DIVIDED DIFFERENCES 


PROBLEMS 

• HARD TO FIND LEADING TERM 

• HARD TO CHANGE INTEGRATORS 

• NO CHEAPER THAN OTHER ESTIMATES 


LOWER ORDER ESTIMATES 

( V t' V tt) 

• ESTIMATE SOLN. GROWTH IN TIME 

• PROS - CHEAP, BETTER THAN GRADIENT ESTIMATES 

• CONS - ACCURATE TRENDS BUT INACCURATE ESTIMATE 

OF MAGNITUDE. 

GRADIENTS 

• USE U x 

PROBLEMS 

• EASY TO FOLL (e.g. FORCING FN. ) 

• NO CHEAPER THAN V t 

• GOOD ONLY FOR SHOCKS 


DEFERRED CORRECTION 


• USES 2 METHODS 

• COMPUTE ERROR ESTIMATE AS A FUNCTION OF THE 2 SOLUTIONS 


PROS 


CONS 


MOST ACCURATE 
ESTIMATES TESTED 


EXTRA TIME LEVELS FOR 2ND METHOD 

DIFFICULT TO FIND 2ND METHOD AND 
ERROR RELATION 


SPECIAL CASE (2h, 

• 2ND METHOD 

• ERROR 


2k) 

USES SAME INTEGRATOR WITH DOUBLE THE STEP SIZES 
V h,k " V 2h,2k 


USE OF DIFFERENTIAL EQ. TO ELIMINATE TIME DERIV. 

• USE U t = f(u,x,t) x TO REPLACE TIME DERIVS. IN TRUNCATION ERROR 
PROBLEMS 

• MESSY TO FIND V m 

• VERY PROBLEM AND METHOD DEPENDENT 

• USEFUL ONLY IF EXTREME PENALTY FOR USING EXTRA TIME LEVELS. 


187 


CONCLUSION 


AUTOMATIC REFINED GR ID GENERATION 

• ARBITRARY ORIENTATION OF RECTANGLES 

• LOW OVERHEAD OF GRID REPRESENTATION 

• REFINEMENTS BASED ON (2h,2k) ESTIMATES 
OF LOCAL TRUNCATION ERROR 


